Nonlinear voltage dependence of shot noise 
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The current noise in a multi-probe mesoscopic conductor 
can have a nonlinear dependence on the strength of driving 
bias voltage. This paper presents a theoretical formulation 
for the nonlinear noise spectra. We pay special attention to 
maintain gauge invariance at the nonlinear level. At small but 
finite voltages, explicit expressions for nonlinear noise spectra, 
expanded order by order in the bias, have been derived. In 
the wideband limit, a closed form solution of the noise spectra 
for finite voltages is obtained. 
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I. INTRODUCTION 

Due to the particle nature of electron, electric cur- 
rent in a conductor fluctuates with time giving rise to 
shot noise. The spectra of shot noise contain informa- 
tion which can be used to characterize electron transport 
in the conductor. For instance it may be used to probe 

and investigate correlations of 
. For this reason shot ngise of 
mesoscopic systems has been studied extensivelytfQ. A 
classical conductor is characterized by Poissonian noiseO, 
where the current fluctuation < (A/)^ > in a frequency 
range Ai/ is proportional to the electrical current /: 



the kinetics of electroE 
electronic wave functio 



< (A/) >= 2q/Az/ where q is the electron charge. In 
a mesoscopic conductor, on the other hand, shot noise is 
also influenced by two other factors: the Pauli exclusion 
principle and the Coulomb interaction. Pauli exclusion 
reduces the classical shot noise by a factor proportional 
to (1 — T) for each transmission subband, assuming the 
transniission coefhcient T to be insensitive to electron 
energytlEI. Coulomb interaction, on the other hand, can 
contribute to reduce or enhance shot noise depending on 
system details. 

The quantum suppression of shot noise has i>^ei|i 
convincingly demonstrated by several experimentfflllj. 
The universal suppression by Coulomb interaction in 
nondegenerate diffusive |Gonductors has been observed 
in computer simulationstil and confixrned theoretically 
using Boltzmann-Langevin equationtj. The quantum 
enhancement of shot noise from the classical value 
due to Coulornb. iateraction has recently been explored 
experimentallyli^O. For a tunneling structure with or 
without a magnetic field, shot noise versus voltage in- 
creases drastically in the region of negative differential 
resistance (NDR). If onfi_assumes sequentialj-tunneling 
of the electron transportllj, numerical resultaij were in 



good agreement with those of the experiment and it in- 
dicated that the enhancement of shot noise was caused by 
Coulomb interaction. From the scattering matrix theory 
point of view, Ref. O examined the effect of Coulomb in- 
teraction and the enhancement of shot noise was found to 
be related to the multistability when a tunneling system 
is out of equilibrium. So far, comparison between exper- 
imental and theoretical resultal3 suggests a need for a 
microscopic theory. 

From a theoretical point of view, focused attention has 
been devoted to cases where the external bias voltage 
strength which drives the current flow is very small. In 
this case the current due to each subband is a linear func- 
tion of the bias voltage: li — {2e'^ /h)VbiasTi, where Ti 
is the zero bias transmission coefhcient of the z-th sub- 
band. Indeed, experiment of Ref. O on quantum point 
contact (QPC) has shown that the peak value of shot 
noise due to first subband is a linear function of Vbtas- 
Similar linear behavior was also found in QPC experi- 
ment of Ref. pG. On the other hand, in many nonlinear 
devices the electric current is a nonlinear function of bias: 
/ = I{Vbias)- A notable example is the resonance tun- 
neling structure where current / varies with bias in a 
nonlinear fashipDj-giszing rise to NDR regions. Previous 
investigationsQoOES indicated that the nonlinear I-V 
characteristics has a profound influence on the shot noise 
spectra, including the enhancement of it. 

The purpose of the present work is to report a micro- 
scopic theory for calculating shot noise at the nonlinear 
regime in mesoscopic conductors. Our theory is based 
on nonequilibrium Green's functions where the electron- 
electron interaction is treated in a self-consistent den- 
sity functional form at the Hartree level. A direct con- 
sequence of self-consistency is that shot noise becomes 
only a function of voltage difference, which is the required 
physical condition (gauge invariance) for a nonlinear the- 
ory. We derive the nonlinear shot noise formula at zero 
temperature, which, in the wideband limit, can be ex- 
actly evaluated. For more general situations we derive 
shot noise spectra order by order in the bias voltage. 

In the next section we present the derivation of the 
nonlinear shot noise spectra. Sections III and IV present 
the wideband limit result with numerical evaluations, as 
well the weakly nonlinear analysis of the shot noise. The 
last section summarizes the main findings of this work. 



II. NONLINEAR SHOT NOISE FORMULA 

The theeretical approachs to shot noiae-jpelude scatter- 
ing matrixQ, semiclassical kinetic-,theoryll3^t2l, nonequihb- 
rium Green's function (NEGF)Ejc3, etc. For a full non- 
linear analysis of quantum transport in the mesoscopic 
regime, we have found that it is most convenient to use 
NEGF, but with the necessary extension of including the 
internal potential build up due to electron-electron inter- 
actions. 

We start from the combined DC thermal and shot noise 
spectra derived in Ref. |^ which is appropriate for a con- 
ductor with non-interacting electrons {Ti = 1), 



< AlaAIfj > 



Av 



q 



xA(l 



J2 / dETr[A^s{a)As^{(3)] 

7(5 

-fs) (1) 



where 



Ai3^{a) = Sa/lSa^ 



^ap'^O'l 



(2) 



where fa — fiE~qVa) is the Fermi distribution function 
and Sap is the scattering matrix. The subscripts denote 
probes which connect our conductor to the reservoirs. 

In order to correctly explore nonlinear voltage depen- 
dence in a quantum transport theory, it is essential to 
include the iking range Coulomb potential as pointed out 
by BiittikeiuJ. It is well known that interacting systems 
are most conveniently dealt within the NEGF formalism. 
For this purpose, we will rewrite Eq.(|I|) using Green's 
functions where the internal Coulomb potential can be 
explicitly included. This way, our nonlinear theory satis- 
fies the gauge invariant condition, where the noise spectra 
remains unchanged when the voltages of all probes are 
shifted by the samepsnstant amount. Recently, Grames- 
pacher and Biittikeiu have discussed the relationship be- 
tween scattering matrix theory and the Hamiltonian ap- 
proach in which the transmission coefficient is expressed 
in terms of Green's function. We will use this approach 
to rewrite Eq. (|l|) in terms of the Green's functions. We 
will then supplement it with the necessary steps of deter- 
mining the internal electro-static potential build-up due 
to Coulomb interactions. 

In the following, we consider a quantum coherent 
multi-probe conductor specified by the following Hamil- 
tonian 

-ff = ^ Ckacl^Cka + Hcen{dn, ^1} 

ka 

+ ^ [Tka^ncl^dn + C.c] (3) 

ka,n 

where £ka = e^ + qVa- The first term of Eq.(g) describes 
the probes where DC signal is applied far from the con- 
ductor; the second term is the general Hamiltonian for 
the scattering region; the last term gives the coupling be- 
tween probes and the scattering region with the coupling 



matrix Tka.n- Here c\^ (cka) is the creation (annihila- 
tion) operator of electrons inside the a-probe. Similarly 
dj, (d„) is the operator for the scattering region. It is im- 
portant to note that we will include the internal Coulomb 
potential U inside the scattering region so that the actual 
Hamiltonian of the scattering region is Hem + qU . 

The retarded scattering Green's function G^ = 
G^{E, U), where U = U{r) is the electro-static potential 
build-up inside the scattering region due to interacting 
electrons, is given bycj 



G''{E,U) = 



(4) 



E-H-qU-T,^ 

where the self-energy S'' = J2a ^ai^ — qVa) is defined 
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and 77 is a positive infinitesimal and Ta{E) is the 
linewidth function 
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The scattering matrix can be expressed. in t erms of 
Green's functions by the Fisher-Lee relationnJc3 



Sa/3 
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where Wq satisfys 2TTWaWl — Fq. Using this relation, 
it is straightforward to show that Eq. (|^) becomesc3, 

< AlaAIp > = ^AiyJ2 f dETriitSasraG"^ 

-iSpsT^G'^^T^G'^TpGnU^il-fs) (8) 

where G° = G°'{E,U) is the advanced Green's function 
and we have used the notation Ta = Ta{E — qVa) for the 
linewidth function. 

Notice that we have explicitly included the internal 
potential landscape U{t) into the Green's functions, but 
this landscape can only be obtained self-consistently. 
This is a crucial step in the development of a gauge in- 
variant nonlinear DC theory. We determine the internal 
potential U{y) by the self-consistent Poisson equation 



V^C/ = Amq f{dE/2TT)G<{E, U) 



(9) 



where the lesser Green's function G^ is related to the 
retarded and advanced Green's functions G'' and G", 

G<{E,U) = G^Y.'Tp{E-qVp)f{E-qVfi)G'' . (10) 

Eq. (p) is, in general, a nonlinear equation because G*"'" 
depends on t/(r) (see Eq.(H)). By self-consistently solv- 
ing Eqs. (p|,p|jlC|) , we obtain the Green's functions as well 



as the internal potential U. Then we can calculate shot 
noise from Eq. (^ which is a general nonlinear function 
of the external bias { Vq}. This theoretical procedure can 
be carried out at least numerically, but in this work we 
are interested in cases where analytical derivations are 
possible. 

For a two-probe system, Eq. (H) reduces to 

< (Alf > = ^Ai.JdE{[hil - h) + /2(1 - /2)] 

xTr[f] + {h^hfTT[{l-f)f]} (11) 

where T{E, U) = riG"'r2G"' is the transmission operator 
such that Tr\T{E, U)] is the transmission coefficient, here 
the trace is over the matrices written in real space. 

To end this section, we discuss the gauge invariance 
condition. It is easy to prove that the noise spectra Eqs. 
(0) and (ni|) are gauge invariant: shifting the potential 
everywhere by a constant Vo^ U -^ U + Vo and Va — > 
Va + Vo, < (A/)^ > calculated from these expressions 
remains the same. It is useful to note that in Eqs. ([^jll|) 
the quantity T depends on bias voltage: without such 
a voltage dependence the gauge invariance can not be 
satisfied. 



III. THE WIDEBAND LIMIT 

In this section we evaluate the shot noise spiectra at the 
wideband limit. In this commonly used limited, the cou- 
pling matrix F is assumed to be independent of energy 
which drastically simplifies the algebra. The wideband 
limit corresponds to cases where the probes have no fea- 
ture, thus the internal potential U{r) becomes a space- 
independent constant Uq (the value of Uq depends on the 
voltages {Va} and it still needs to be determined). For a 
single level system, as far as the nonlinear current-voltage 
curve is concerned, the wideband limit corresponds to a 
resonance tunneling system where the scattering matrix 
has the Breit-Wigner form. 

In wideband linfit the steady state Green's function 
takes a very simple form, Gg = 1/{E — Eq + iT/2), thus 
the integral in Eq. (O) can be done exactly at zero 
temperature. We obtain. 
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where AEf3 = Ep - Eq - qUo + qVp . 

While the internal potential Uq can be determined by 
the Poisson equation (0), which requires numerical anal- 
ysis, we instead determine it by introducing the geomet- 
rical capacitances Ci and C2 of the left and the right 



coupling regions (regions where our conductor connects 
to the two probes) respectively. The charse, in the well 
due to the Coulomb interaction is given byE3 



AQ = -I / {dE/2n)[G<{E,UQ) ~ 
^Ci{Uo^Vi) + C2{Uo-V2) 



Go<] 



(13) 



where AQ is the total charge in the well and Gq is the 



equilibrium lesser Green's function, 
limit, this equation reduces to 



In the wideband 
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When Ci = C2 = 0, Eq.([l4|) corresponds to the quasi- 
neutrality approximation which neglects the charge po- 
larization in the system. Thus using two phenomenolog- 
ical constants Ci and C2 we can determine Uq from the 
last equation. Hence the noise spectra of Eq. ( p^ ) is now 
completely specified. Finally, we can check the gauge 
invariance of Eq. (O) : the wideband approximation and 
the charging model for Uq do not disrupt the satisfaction 
of it. Indeed, raising both V3 and Uq by the same amount 
does not alter the shot noise given by Eq. (|l2). 

In Fig.l, we have plotted the differential shot noise 
d < (A/)^ > /dV as a function of the voltage for 
four different set of system parameters: symmetric struc- 
tures with Fi = F2 = 0.5, Ci = C2 = 0.5 (solid 
line); Fi = F2 = 0.5, Ci = C2 = 0.1 (dotted line); 
Fi = F2 = 0.1, Ci = C2 = 0.4 (dot-dashed line); and an 
asymmetric structure with Fi = 0.1, F2 = 0.8, Ci = 0.1, 
and C2 = 0.8 (dashed fine). For symmetric structures, we 
always observe two peaks for the differential shot noise 
whereas for the asymmetric structure there is only one. 
This can be understood qualitatively as follows. Since the 
shot noise at small voltage is proportional tcO T — T^, the 
suppression reaches maximum near the resonant point 
for the symmetric structure because T w T^ « 1 at reso- 
nance. As a result, a peak appears on each side of the res- 
onant point giving rise to the two peaks in Fig. 1. For an 
asymmetric structure, however, the shot noise suppres- 
sion is not as strong as that of the symmetric case (see 
Fig. 2). For the very asymmetric structure used here, the 
quantum resonance is very weak, resulting to only one 
peak in the differential shot noise spectra. For the sym- 
metric cases, the separation between two peaks is pro- 
portional to F. The different resonant positions for solid 
line and dotted line in Fig. 1 (the dips near voltage equals 
4 and 6 respectively) is due to the Coulomb interaction 
in the well. Smaller capacitance coefficients (dotted line) 
correspond to large internal potential, which shifts the 
level position from Eq to higher values Eq + qUq. This 
is why that the resonant position for the dotted curve 
(with smaller capacitance coefficients) is shifted further 
relative to the solid curve (with larger capacitances). 



As discussed in the Introduction, the classical shot 
noise is given by 2qIlS.v. The deviation from this classical 
value is usually characterized by the Fano factor which 
is defined as 



7 



< (A/)' > 
2qIAiy 



(15) 



In the wideband limit, it is easy to derive the current / 
to be, 
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Before presenting the plot of 7 for our nonlinear analysis, 
two observations are in order. First of all, it is easy to 
show that in the limit {Vi — V2) -^ and F — » 0, the 
Fano factor 7 — > 1. Secondly, in the opposite limit when 
voltage difference is very large_the Fano factor is given 
by the well known expressior£3 7 — (Ff + F2)/F^. For 
symmetrical systems (Fi = F2) this large voltage limit 
of 7 is 0.5|— I The same behavior has been observed in 
experimentll3 except near the NDR region. Physically, 
the currents from both leads contribute to the Fano fac- 
tor. At large voltage, the current from the low biased 
lead can be neglected and the shot noise is suppressed. 

In Fig. 2, we have depicted the Fano factor versus 
voltage. The system parameters are the same as that of 
Fig. 1. As expected, the Fano factors approach to 0.5 for 
symmetric structures at large voltage. We also see that 
the Fano factors are minimum near the resonance. In 
the wideband limit, the Fano factor can be smaller than 
0.5 for the symmetric case. For smaller F, the transition 
from 7 « 1 to 7 ~ 0.5 is much sharper (dot-dashed line). 
Hence this result suggests a more pronounced noise re- 
duction for conductors which are more weakly coupled to 
the leads (smaller F). On the other hand, for the asym- 
metric case the suppression of the Fano factor is not as 
strong (dashed line). This is due to the fact that quan- 
tum resonance is not as well established in an asymmetric 
system as that in a symmetric one. 

The reason that our full nonlinear results of Fig. 2 
only shows shot noise reduction in a resonance system is 
due to the fact that we have applied the wideband limit 
for the Green's function. Since wideband limit does not 
allow negative differential resistanceo, we can not ob- 
serve NDR and hence the enhancement of shot noiseli^. 
To obtain NDR within the wideband approximation, we 
assume the leads to have a finite occupied bandwidth by 
introducing an energy cutoff itLthc integration of Eq. ( 12 ) , 
as suggested by Jauho et a£3. A result of this simple 
procedure indeed produced a Fano factor which can be 
greater than unity, as shown in Fig. 3. Our analysis thus 
reconfirms that shot noise can be enhanced by the ex- 
istence-of a NDR region. However, the experimental 
resultal3 showed a much sharper increase of the Fano 
factor when bias is varied than that showed by Fig. 3. 
It thus seems that one needs to go beyond the wideband 
limit to obtain a detailed quantitative agreement. 



IV. WEAKLY NONLINEAR LIMIT 

The wideband limit discussed above reduces the sys- 
tem to essentially a single level and zero-dimensional 
quantum dot, but this allows us to obtain closed form 
results for the full nonlinear shot noise spectra including 
large bias voltages. In this section we examine the an- 
other limit, namely the weakly nonlinear limit where the 
bias is finite but not large. In this case we can expand 
the shot noise formula order by order in bias, and derive 
the weakly nonlinear shot noise spectral coefficients. 

For small bias voltages, we expand the noise spectra of 
two-probe systems Eq. (O) in terms of it. 



< (A/)2 >= Pa + P,{Vi - V2) + P2{Vi 



V2? 



(16) 



where the equilibrium noise Pq and the linearjioise spec- 
tra P\ have been considered in detail beforeQ. Here we 
derive an expression for the second order non-linear noise 
spectra P2 at zero temperature. To proceed, we first 
need to determine the internal potential [/. In the last 
section we applied a phenomenological nonlinear capac- 
itance charging model to find C/, here we will solve V 
self-consistently order by order in the bias. We expand 
the internal potential V in powers of voltages. 



C/ = C/e, 



L^ 
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where Ueq is the equilibrium poteptiaLand Ua{v), Ua^..(r) 
are the characteristic potentialsE3E2l. They are deter- 
mined by the Poisson like equations which are obtained 
by expanding Eq. (0) in powers of voltage, 

.dna{x) 
dE 
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Il{x,x') is the Lindhard functionE3'E2l defined as 

U{x,x') = -i l{dE/2TT)f[Gl{x,x')Gl{x',x) 



Gl{x,x')G''a{x\x)] 
'^no(x) 
q5U{x') 
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where / = f{E), Gq is the equiUbrium Green's function, 
and the generating function 1X0(2;) is given by 

Uo{x) = ~i j{dE/2^)f[Gl{x, x) - Gg(x, x)] (21) 

and 5/5U{x2l is the functional derivative defined in Ref. 
|33| . In Eq.([l9|), Ii{x,x' ^x") is the second order nonhnear 
response function defined as 



n(x,x',x" 
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and n„ the Lindhard function for lead a 
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II^{x,x') = / (d^/2^)/[(GSr„Gg),,-(Gg),-, + CO.- 



(23) 

so that n(x^') = J2a^ot{x,x'). The partial local den- 
sity of stataHd at contact a, dua/dE, called the injectivity 
is given by, 



dna{x)/dE = Ila{x,x')dx' 



(24) 



and d^Ua/dE'^ is the energy derivative of the injectiv- 
ity. Finally, dn/dE = J2a^'^a/dE is the local DOS. It 
can be shown that tJaCjCharacteristic potential satisfy the 
following sum rulesEETEHl, 



^Ua^l 



and 
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4/3}i 



(25) 



(26) 



Here the subscript {/3}/ is a short notation of I indices 
7,(5,77,- ■ •. 

With these prepartions we can now derive the second 
order nonlinear shot noise coefhcient ^2- At zero tem- 
perature, the shot noise formula Eq. (pi]) reduces to 

2 pEp+qVi 

< (A/)2 >= Aiy— / dETr[{l - f)f] . (27) 

Denoting g{E, U) = {I - f)f, we expand g{E, U) with 
respect to E and V: 
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where U is the diagonal matrix for the internal potential, 
50 = 9{Ef,0) and idgo/dU)U = j:^(S^/SU{x))U{x). 
Substitute the above equation into Eq. (p7|) and complete 



the integral over energy E, the noise spectra up to the 
second order in voltage is 
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Using the relationtil qdgo/dE = —dga/dU for gauge in- 
variance and Eq.(p3), we arrive at 



''-^''li^-l|7(2«--i)l 



(30) 



Following the same line of development, we can de- 
rive higher order nonlinear shot noise coefficients. For 
instance, the third order nonlinear noise spectra is found 
to be, 

P, = Ai.fTr[3q^u,, + ^(1 - 3u, + 3ul)] . (31) 
OTT dU dU'^ 

As an explicit example, let's derive P2 and P3 for a res- 
onance tunneling system using scattering approach and 
comparing 'vrith the result of NEGF. Using the Breit- 
Wigner formcil for the scattering matrix near a resonance 
energy E^, Sap{E) ~ [Sap - i^FaF^/A], where F^ is 
the decay width of barrier a, A = E — Eq + iT/2 with 
F = Fi + F2, we obtain. 
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and 
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where T — F1F2/IAP. In the derivation of Eqs.(|32)-|and 
(p3|), we have used the quasi-neutrality conditioiE^ for 
determining the characteristic potential so that ui — 
Fi/F and un == -2{E ~ Eq)T/T^. Since resonance 
tunneling with Breit-Wigner scattering matrix is equiv- 
alent to the wideband limit of the last section, expres- 
sions (K2 33 ) can be directly obtained by expanding the 
wideband limit results ( |l2| , [l4|) to the appropriate order in 
voltage. It is straightforward to prove that the same re- 
sults are obtained from this direct expansion. This gives 
a confirmation on the validity of the weakly nonlinear 
analysis presented here. 



V. SUMMARY 

In this work, we have developed a general nonlinear DC 
theory for calculating the shot noise spectra in the meso- 
scopic regime. The framework is based on on nonequilib- 
rium Green's functions with the important extension of 



solving the internal potential build up self-consistently. 
A direct advantage of our method is that the final ex- 
pression for shot noise becomes gauge invariant which 
is an essential requirement for any nonlinear transport 
theory. Eqs. dq,0,@) completely determine the nonlinear 
shot noise spectra of an arbitrary multi-probe conduc- 
tor, they form the basic results of our theory. Practi- 
cally, one must solve the quantum scattering problem 
which gives the Green's functions, in conjunction with 
the Poisson equation. Technically these expressions form 
a convenient basis for numerical predictions of shot noise 
spectra at finite bias voltages. For instance one can eas- 
ily compute various Green's functions and the coupling 
matrixptj for multi-probe conductors using tight-binding 
modelscJ; and the Poisson equation can be solved in real 
space using very powerful numerical techniqueg23. 

In the wideband limit and the weakly nonlinear limit, 
the basic equations (HJj,^) can be analyzed in closed 
form. Our nonlinear theory reveals that the shot noise 
of a mesoscopic conductor can be quite sensitive to the 
external bias strength, and in general the suppression of 
noise is most efficient near a quantum resonance point, 
and is stronger for symmetric systems than asymmetric 
ones. The suppression is also more efficient for conduc- 
tors weakly coupled to the leads. In the presence of nega- 
tive differential resistance region of the nonlinear current- 
voltage characteristics, our result confirms the existence 
of shot noise enhancement which has been observed ex- 
perimentally. For weakly nonlinear transport regime, we 
have derived the shot noise nonlinear coefhcients order by 
order in bias, and these coefhcients should be adequate 
when the external bias is finite but not large. 
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FIGURE CAPTIONS 

Fig. (1) The differential shot noise versus the voltage. Solid 
hnc: Fi = Fa = 0.5 and Ci = C2 = 0.5; dotted line 
Fi = F2 = 0.5 and Ci = C2 = 0.1; dot-dashed line 
Fi = F2 = 0.1 and Ci = C2 = 0.4; and dashed line 



Ti = 0.1, Ta = 0.8, Ci = 0.1, and C^ == 0.8. Here 
Ef- Ea^ -2.0. 

Fig. (2) The corresponding Fano factor of Fig. (1). 

Fig. (3) The Fano factor versus the voltage when the energy 
cut off is introduced for Fi = r2 = 0.5, Ci = C2 = 
0.5 and Ep - Eq = -2.0. 
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